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Abstract 

The scattering of a straight, infinitely long string by a rotating black hole is 
considered. We assume that a string is moving with velocity v and that initially 
the string is parallel to the axis of rotation of the black hole. We demonstrate 
that as a result of scattering, the string is displaced in the direction perpendic- 
ular to the velocity by an amount k{v, b), where b is the impact parameter. The 
late-time solution is represented by a kink and anti-kink, propagating in oppo- 
site directions at the speed of light, and leaving behind them the string in a 
new "phase" . We present the results of the numerical study of the string scat- 
tering and their comparison with the weak-field approximation, valid where 
the impact parameter is large, b/M ^ 1, and also with the scattering by a 
non-rotating black hole which was studied in earlier works. 
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1 Introduction 



Study of cosmic strings and other topological defects and their motion in an external 
gravitational field is an interesting problem. Cosmic strings are topologically stable 
one-dimensional objects which are predicted by unified theories. Cosmic strings (as 
well as other topological defects) may appear during a phase transition in the early 
Universe. A detailed discussion of cosmic strings and other topological defects can be 
found in the book by Shellard and Vilenkin Cosmic strings are naturally predicted 
by many realistic models of particle physics. The formation of cosmic strings helps 
to successfully exit the inflationary era in a number of the inflation models motivated 
by particle physics ^. The formation of cosmic strings is also predicted in the 
most classes of superstring compactification involving the spontaneous breaking of a 
pseudo-anomalous U{1) gauge symmetry ( see e.g. 0] and references therein). Recent 
measurements of CMB anisotropy and especially the position of the acoustic peaks 
exclude some of the earlier proposed scenarios where the cosmic strings are the main 
origin of CMB fluctuations. On the other hand, recent analysis shows that a mixture 
of inflation and topological defects is consistent with current CMB data ||, 0, H, . 

A black hole interacting with a cosmic string is a quite rear example of interaction 
of two relativistic non-local gravitating systems which allows rather complete analysis. 
This makes this system interesting from pure theoretical point of view. From more 
"pragmatic" viewpoint, this system might be a strong source of gravitational waves. 
But in order to be able to study this effect one needs first to obtain an information 
on the motion of the string in the black hole background. The situation here is very 
similar to the case of gravitational radiation from bodies falling into the black hole. 

In this paper we consider the scattering of a long cosmic string by a rotating 
black hole. We assume that a string is thin and its gravitational back reaction can 
be neglected. If r/ is a characteristic energy scale of a phase transition responsible for 
a string formation then the thickness of a string is p ~ r^"^ while its dimensionless 
mass per unit length parameter /i* = Gfi/c^ ~ rj^. For example, for GUT strings 
Pgut ~ 10~^^cm and p^fyT ~ 10~^, and for the electroweak phase transitions pew ~ 
10~^^cm and /x^^ ~ 10~^^. Since /x* -C 1 one can neglect (at least in the lowest 
order approximation) effects connected with gravitational waves radiation during the 
scattering of the string by a black hole. 

We assume that the size of the string is much larger then the Schwarzschild 
gravitational radius rs = 2GM/c^ of a black hole and its total mass is much smaller 
than the black hole mass. The latter condition together with /i* <C 1 means that we 
can consider a string as a test object. In order to specify the scattering problem we 
consider the simplest setup, with the string initially far away from the black hole, 
so that the string has the form of a straight line. We assume that string is initially 
moving with the velocity v towards the black hole and lying in the plane parallel to 
the rotation axis of the black hole. This plane is located at distance h with respect 
to the parallel "plane" passing through the rotation axis. In analogy with a particle 
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scattering we call b the impact parameter. 

The scattering of a cosmic string by a non-rotating black hole was previously 



studied both numerically and analytically [10, 11, O, O, nM. For b ^ rs the string 



moves in the region where the gravitational potential GM/r is always small, and one 
can use a weak-field approximation where the string equation of motion can be solved 
analytically Jll], |12|, . 

String scattering in the weak field limit can be described in qualitative terms. 
In the reference frame of the string the gravitational field of the black hole is time 
dependent and it excites the string's transversal degrees of freedom. This effect occurs 
mainly when the central part of the string passes close to the black hole. Since 
information is propagating along the string at the velocity of light, there always exist 
two distant regions, 'right' and 'left', which have not yet felt this excitation. This 
asymptotic regions of the string continue their motion in the initial plane ('old phase'). 
After scattering, when string is moving again far from the black hole, its central part 
moves again in a plane which is parallel to the initial plane but is shifted below it 
by a distance k = 2110 Mv /y/c^ — v"^ ('new phase'). There exist two symmetric kink- 
like regions separating the 'new' phase from the 'old' one. These kinks move out 
of the central region with the velocity of light and preserve their form. Besides the 
amplitude k, the kinks are characterized by a width w which depends on the impact 
parameter and the initial velocity of the string. 

In the strong-field limit this picture qualitatively remains the same for string scat- 
tering by a Schwarzschild black hole. This problem was analyzed in detail numerically 
in |[Tl|. In particular, the quantities k, and w were found for different values of b and 
V. The numerical calculations in the strong field limit demonstrated also that for 
special values of the initial impact parameter b < bcritiv) a cosmic string is captured 
by the black hole. The critical impact parameter for capture by the Schwarzschild 
black hole was obtained in ||1^, |13|. 



In this paper we study of scattering of straight cosmic strings by a rotating black 
hole. We assume that the initial direction of the string is parallel to the rotational axis 
of the black hole. Our focus is on the calculation of the same quantities, k and w, and 
to study how black hole rotation modifies the earlier results for a non-rotating black 
hole. We present here results for string spacetime evolution, real profiles of a string 
at different external observer times T, as well as asymptotic scattering data. Results 
connected with capture of the string and near-critical scattering will be discussed in 
another paper. 

The paper is organized as follows. Equations for string motion in a curved space- 
time and their solutions in the weak field regime are discussed in section 2. Section 
3 contains the formulation of the initial and boundary conditions for straight string 
motion in the Kerr geometry, and the general scheme of the numerical calculations. In 
section 4 we describe results for strong field scattering of strings in the Kerr spacetime 
and "real time" profiles of the strings. Section 5 contains information about asymp- 
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totic scattering data (displacement parameter and kink profiles) and their dependence 
on the impact parameter, velocity, and angular velocity of the black hole. Section 6 
contains a general discussion. Numerical details are supplied in an appendix. 



2 Equations of motion 
2.1 String equation of motion 

The worldsheet swept by the string in (4D curved background) space-time can be 
parametrized by a pair of variables {A = 0, 1) so that string motion is described 
by the equation = X^{(^). The dynamics of the string in a spacetime with metric 
g^i, is described by the Nambu-Goto effective action 

= -^Jd'C^, (2.1) 

where 

is the induced metric on the worldsheet. The same equations of motion of the string 
can be derived from the Polyakov's form of the action |15|, 



hAB] = / d'cV^h^'^^AB . (2.3) 



We use units in which G = c = 1, and the sign conventions of [jT6|. In ( |2.3| ) Hab is 
the internal metric with determinant h. 

The variation of the action ( |2.3| ) with respect to and Hab gives the following 
equations of motion: 

nX" + h'^^T^fsX^^Xj^j^ = , (2.4) 
7AB-lhABh^'''7cD = 0, (2.5) 



where 



□ = -^^A{^/^h^''^B) . (2.6) 



The first of these equations is the dynamical equation for string motion, while the 
second one plays the role of a constraint. 

We fix internal coordinate freedom by using the gauge in which Hab is conformal 
to the flat two-dimensional metric tjab = diag(— 1, 1). In this gauge the equations of 
motion for the string have the form 

nX^" + r^'ap X'X X% 7]^^ = , (2.7) 
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and the constraint equations are 



701 = 9^v-q:^-q^ = , (2.8) 

fdX^'dX'' dXf'dX''\ , , 

7«o + 1.>=9..(^^ + ^^j =0. (2.9) 

Here, t = ^ o = C} and □ = —9^ + 9^. These constraints are to be satisfied for 
the initial data. As a consequence of the dynamical equations they are valid for any 
later time. 



2.2 Weak field approximation 

In the absence of the external gravitational field g^^y = rjf^^, where 77^,^ is the fiat 
spacetime metric. In Cartesian coordinates (T, X,Y, Z), rj^i, = diag(— 1, 1, 1, 1) and 
= 0, and it is easy to verify that 

= X^{T,a) = (rcosh/5,rsinh/? + Xo,6,a) , (2.10) 

hAB = Vab = diag(-l, 1) , (2.11) 

is a solution of equations (|2.4|) and (|2.5| ). This solution describes a straight string 
oriented along the Z-axis which moves in the X-direction with constant velocity 
V = tanh/?. Initially, at tq = 0, the string is found at X^{0,a) = (0,Xo,6, cr). Later 
when we use this solution as initial data for a string scattering by a black hole, b will 
play the role of impact parameter. For definiteness we choose 6 > and Xq < 0, so 
that /3>0. 

Let us consider how this solution is modified when the straight string is moving 
in a weak gravitational field. We assume 

Qixv = 'nf,u + q^,y , (2.12) 

'^'^(C)='^'o"(C) + X^(C), (2.13) 

where is the metric perturbation and is the string perturbation. By making 
the perturbation of the equation of motion one gets 

Ux^ = r, (2.14) 

where 

r = no = -r^f^iXo) x,yx,%v^^ . (2.15) 

We use Cartesian coordinates here so that F^q^ is simply the Christoffel symbol for 
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The linearized constraint equations are 



As in the exact non-hnear case, if these hnearized constraints are satisfied at the initial 
moment of time r they are also valid for any r for a solution of the dynamical 
equations ( p.l4| ). 

The linearized equations can be used to study string motion in the case where it is 
far away from the black hole. In this case the gravitational field can be approximated 
as follows: 

ds" = -{l-—)dT^ + {l + —){dX^ + dY^ + dZ^)-%{XdY-YdX)dT, (2.18) 
R R R'^ 

where R^ = X"^ + + Z"^, and M and J = aM are the mass and angular momentum 
of the black hole, respectively. In fact this asymptotic form of the metric is valid for 
any arbitrary stationary localized distribution of matter, provided that the observer 



is located far from it. In agreement with (|2.18| ) the gravitational field perturbation 
can be presented as 

q^.u = qj!u + qil, (2.19) 
where the Newtonian and Lense-Thirring parts are 

g;X = 2^5^,, ^=^, g^J = e,),03X" . (2.20) 

Here eajj-ys is the antisymmetric symbol. The Lense-Thirring part g^J of the metric 
T7| is produced by the rotation of the source of the gravitational field and it is 



responsible for frame dragging. 



2.3 Newtonian scattering 

In the linear approximation we can study the action of each of the parts of the 

metric perturbations independently. For the Newtonian part one has the following 
expression for the force 

2 sinh(/5) cosh(/5) , (^,1 (2.21) 

0, (2.22) 

-2 sinh^ /3(^_2 , (2.23) 

-2 cosh^ /3(^,3 , (2.24) 



J N — 

In = 

J N — 

In = 
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and the constraint equations read 

xl'.-cosh/^x^. + sinh/Jx;. = 0, (2.25) 
x\ - cosh /?x°r + sinh (3x]r = -'^f cosh^ /3 . (2.26) 

The result = means that in the first order approximation the string hes in 
the yZ— plane at any fixed time r = const. The dynamical equations ( p.3|) can be 
easily integrated by using the retarded Green function (see ||TT[] ). 



We focus now our attention on the behavior of which describes the string 
defiection in the plane orthogonal to its motion. It is easy to show that the asymptotic 
value of x^(t ~^ oo, a) is the same for any fixed a. We denote it b — k. The following 
expression is valid for k 

dXdZf, (2.27) 



2 sinh/? 

where Ho is the two-dimensional worldsheet swept by the string in its motion in the 
background spacetime. This integral can be easily calculated since 

dXdZ^y = 2nM (2.28) 

no 

is the fiux of the Newtonian field strength through Ho. Thus we have 

K = 27rM sinh /5 . (2.29) 

At any late moment of time r the central part of the string is a straight line moving 
in the plane Y = b — k, while its far distant parts move in the original plane Y = b. 
These different "phases" are connected by kinks propagating in the direction from the 
center with the velocity of light. (For details, see |13|). We call k the displacement 
parameter. Figure |l] shows the Y-direction displacement for the Newtonian scattering 
as a function of r and a. 

2.4 Lense-Thirring scattering 

In the presence of rotation the Lense-Thirring force acting on the string in the lin- 
earized approximation is 

/o, . 6Jsinh^/3M^^^^, (2.30) 
/It = 0, (2.31) 
fl^ = 2Jsinh/?cosh/3(^^-|^^ , (2.32) 

flj. = 6J sinh/? cosh (2.33) 
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Figure 1: F— profile for Newtonian scattering 



wfiere 7^ = ^(XoTTshlh^Sj^TFTo^. As before, f = 0. 
The constraint equations now take the form 

- cosh fJxl + sinh Px'. = , (2.34) 

X^^ - cosh /3x°r + sinh Px\ = -2 J sinh P cosh . (2.35) 

The dynamical equations can be solved analytically. For the initial conditions 

X°|r=o = 0, (2.36) 

X'|r=0 = ^|r=0 = 0, (2.37) 



dr 

X'|r=0=^|r=0 = 0, (2.38) 



the solution is 



n T7 1 /I / (-^0 + ''■sinh /3) sinh /3 + cr (Xq + r sinh /?) sinh /3 — o" 
X = Jbsmh(3 h 



(+-)'R- ^(++)^ 
sinh (3 — {t — a) Xq sinh /9 — (r + cr) 



, > (2-39) 
X' = 0, (2.40) 
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P , . , ^ , ^ f t"^ sinh^ 13 + Ta sinh^ (3 + 2Xor sinh 3 + cxXq sinh /3 + Xn^ + 6^ 
Y = J smh o cosh o — 

-r^ sinh^ 13 + ra sinh^ /5 - 2Xor sinh {3 + aXo sinh [3 - Xl-lP' 
^ Xo sinh /5(r + a) + + 52 X^ sinh /5(r - a) + X^ 



^{++)P+ ^{—)P- 

., TT • , /I , „/ (Xo + r sinh /5) sinh /3 — 0" (Xq + r sinh /9) sinh /3 + a 
X = Jb smh p cosh p | — — 



V A(++)7^ A(__)7^ 



Xo sinh (3-{T + a) ^ Xpsinh/? - (r - d) 



^±± = (^0 ± (r ± a) sinh + 6^ cosh^ p , (2.43) 



P± = V^o + + (r ± a)2 . (2.44) 

Using these relations it is possible to see that the displacement parameter n = 
lim^^oo x'^i'^j'^ = 0) vanishes for this solution. The reason is the following. Using 
( |2.32| ) it is easy to present flj, as 

fl^ = 2 J sinh /? cosh/5— (^) . (2.45) 

The asymptotic displacement is given by ( |2.27| ). The integral over Z or, equivalently, 
the integral over a from which is a total derivative over a, reduces to boundary 
terms a/TZ^ at a = ±00, which vanish. 

Figure ^ shows the y— direction displacement as a function of r and a for the 
weak field scattering. It should be emphasized that the scale of structures for Lense- 
Thirring scattering is much smaller than that for Newtonian scattering. This can be 
easily seen if we compare the Newtonian force Jn ~ M/TZ'^ with the Lense-Thirring 
force J'lt ~ J/TZ^ 

ll^L^J_<J_, (2.46) 
fN Mil - Mb ^ ' 

For the scattering by a rotating black hole J = aM, where \a\/M < 1 is the rotation 
parameter. Hence 

frr M 

^ < - (2.47) 
In 

which is small for the weak field scattering. For this reason, in the weak field regime 
the string profiles for the prograde or retrograde scattering by rotating black hole 
do not greatly differ from the profiles for the scattering by a non-rotating black 
hole of the same mass. The situation is different for strong-field scattering: the 
displacement parameter is different for prograde and retrograde scattering, and for 
near-critical scattering the profiles of the kinks contain visible structure produced by 
effects connected with the rotation of the black hole. 
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Figure 2: F— profile for Lense-Thirring scattering 



3 String scattering by a Kerr black hole 
3.1 String in the Kerr spacetime 

Our aim is to study a string motion in the Kerr spacetime. The Kerr metric in 
Boyer-Lindquist coordinates {t, r, 6, 0) has the form 

ds^ = -(1 )df + —dr^ + J:de^ + rf^^ dtd(f) , (3.1) 

S = + (3.2) 
A = - 2Mr + , (3.3) 
A = {r^ + a^f -a^Asin^e , (3.4) 

where M is the mass of the black hole, and J = aM is its angular momentum 
(0 < a < M). At far distances this metric has the asymptotic form (|2.18| ). 

In order to be able to deal with the case where part of the string crosses the 
event horizon for the numerical simulation we adopt the so-called Kerr (in-going) 
coordinates {v, r, 9, cf)) 

,2 2Mr, 4aMrsin2^,^ ~ .2^, ,7 
ds = — (1 :^)dv + zdvdr — dvd(p — asm Udrdcp 

^ ,^0 y4sin^ , ^ 

+ T.d6^ + — - — d<p^ , (3.5) 
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with 



t + r + Mln 



+ 



In 



+ M 



2VM2 - a2 



In 



- M - VM2 - a2 



- M + - a2 



(3.6) 
(3.7) 



The metric ( p.l| ) has the asymptotic form 



(is^ = -dt^ + dr^ + r\d9^ + sin^ 



+ 



2M 



AaM 



sm^edtdcj). (3.8) 



Let us introduce new "quasi- Cartesian" coordinates 



T 
X 
Y 
Z 



where 



R sin 6^ cos ( 
R sin 6' sin <; 
Rcos9 , 



R + M. (3.9) 
in the quasi-Cartesian coordinates (3.9) 



One can easily check that the metric 
has the asymptotic form (|2.18|) 

It should be emphasized that without the shift (|3.9| ) of the radial coordinate one 
does not recover ( |2.18| ). The reason for this is easy to understand if one considers the 
same problem for the Schwarzschild geometry. The asymptotic limit of (|2.18|) with 
a = can be found in the isotropic coordinates with r = R{1 + M/2R)^ in which 
the spatial part of the metric is conformal to the flat metric. Asymptotically, it is 
sufficient to use ( |3.9D which is the leading part of r at large R. 

In what follows we shall use these quasi-Cartesian coordinates for representing 
the position and the form of the string in the Kerr spacetime even if we are not 
working in the weak field regime. It should be emphasized that since the space in the 
Kerr geometry is not fiat, plots constructed in these coordinates do not give a "real 
picture" . This is a special case of a general problem of the visualization of physics in 
a 4-dimensional curved spacetime. 



3.2 Initial and boundary conditions 

In an ideal situation, we would study scattering by initially placing our infinitely 
long string at spatial infinity, where the spacetime is fiat and the straight string can 
be described in simple terms. In a numerical scheme, however, the string cannot be 
infinitely long and we must start the simulation at a finite distance from the black 
hole. We discuss here the initial and boundary conditions used for the simulations. 
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In studying the scattering of a straight string, we consider the special case where 
the string is initially parallel to the axis of rotation of the black hole. We use Cartesian 
coordinates in the asymptotic region so that X-axis coincides with the direction of 
motion and the Z-axis is along the string (see Fig. ^). 



We consider a string segment of length L ~ lO^M at a time Tq and an initial 
distance Xq from the black hole. In order to keep accuracy high and yet prevent the 
calculation time from being inordinately large, we do not evolve the straight string 
numerically from this initial position. Instead, we use the weak-field approximation 
to describe the string configuration at a later time where the distance Xs is closer 
to the black hole, -C |Xo| Q. 

Given a sufficiently long string segment, the boundary points move at a great 
distance from the black hole, so that their evolution can be described by the weak field 
approximation until information about the interaction with the black hole reaches 
the boundary. We denote this time as r*. Starting from this moment we solve the 
dynamical equations in the region \a\ < L — {t — t^) (see figure The larger is the 
initial value L the longer one can go in r in the simulation. We choose L to be large 
enough to provide the required accuracy in the determining the final scattering data. 

Since the boundary conditions found from the weak field regime calculations are 
not exact there will be disturbance created by this effect. In some cases, when it is 
important to exclude these disturbances, we used a modified scheme of calculations 
in which no boundary conditions are used from the very beginning of the simulation. 
The price for this is that one needs to take longer initial size of L which in turn 
increases the computational time. 

^ In our simulations we put Xq — — KfM and — — 500A/ 
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Figure 3: Initial setup of the scattering experiment 
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END OF SIMULATION 
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Figure 4: Scheme of time domains for scattering problem 



3.3 Solving dynamical equations and constraints 

Using the initial conditions at and boundary conditions we solve numerically the 
dynamical equations (pTTI). Since the equations are of second order, we use the weak- 
field approximation to get initial data at and + At in order to completely 
specify the initial value problem. The numerical scheme uses second-order finite 
differences and evolves the string configuration using an implicit scheme. Because of 
the symmetry a —a it is sufficient to evolve only half of the string worldsheet, 
cr < 0, and use a reflecting boundary condition at the string midpoint. The numerical 
grid is non-uniform with the denser part following the motion of the kink. We used 
the constraint equations (|2.8[ - p.9|) for an independent check of the accuracy of the 
calculations. More details on the numerical scheme can be found in the Appendix. 
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Weak field b=80M scattering (a/M=0) 



Strong field b=6M scattering (a/M=0) 
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4 String profiles for strong field scattering 



4.1 General picture 

Figure |^ demonstrates general features of straight string scattering, here for a string 
with velocity v = 0.762 [jS = 1). The four panels show the displacement in the 
y— direction as a function of internal coordinates (r, a) for weak and strong field 
scattering. 

Figure |^ (a) shows the scattering for an impact parameter 6 = 80 M. The value 
of the displacement parameter at the largest r shown is k = 7.23M. This value 
as well as the form of the kinks are in a very good agreement with the weak-field 
scattering results. Plots for prograde and retrograde scattering at 6 = 80 M are 
virtually indistinguishable from Fig. ^ (a). 

Figures | (b-d) show the scattering for b = 6M. Qualitatively these plots are 
similar to weak-field scattering. The main differences are the following: (1) The 
value of the displacement parameter k = 11.03M for the non-rotating black hole is 
significantly larger then k = 7.23 for b = 80 M. (2) Effects of rotation are more 
pronounced. The displacement k, for retrograde scattering is essentially larger than 
K for the prograde scattering. (3) The string profile after scattering is more sharp, 
the width of the kinks is visibly smaller that for weak-field scattering. 

4.2 "Real-time profiles" of the string for strong-field scat- 
tering 

It should be emphasized that figure |^ give an accurate impression of the displacement 
effect, but they give a distorted view of the real form of the string. The reason is 
evident. The grid imposed on the worldsheet is determined by the choice of (r, cr) 
coordinates. But a r = const section differs from a time T = const section in the 
laboratory slice of spacetime. The position of the string at given time T can be found 
by using the functions T(r, a), X(r, cr), Y{T,a), Z{T,a), to find functions X(T,Z) 
and Y{T, Z). For fixed T these functions determine a position of the string line in 
3-space. 

Figures ^ and || show a sequence of "real-time" profiles for strong-field string 
scattering at five different coordinate times T. We ordered them so that the larger 
number labeling the string corresponds to later time. In figure |^ the black hole is non- 
rotating. Figures |^ and |^ show the "real-time" profiles for retrograde and prograde 
scattering for a maximally rotating black hole, respectively. Again we can see that 
due to frame dragging, the distortion of the string is more pronounced for retrograde 
scattering. 
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Figure 6: "Real-time profiles" of the string and their XZ— and Y Z— projection for 
the strong field scattering by the Schwarzschild black hole. 
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Figure 7: "Real-time profiles" of the string and their XZ— and Y Z— projection for 
the strong field retrograde scattering by the extremal Kerr black hole. 
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Figure 8: "Real-time profiles" of the string and their XZ— and Y Z— projection for 
the strong field prograde scattering by the extremal Kerr black hole. 

5 Late time scattering data 
5.1 Displacement parameter 

Figure ^ shows the dependence of the displacement parameter k on the impact pa- 
rameter h for different velocities v. For a given value of h the curve for a/M = — 1 
always lies higher than the curves for more positive values of ajM. This is true for 
all values of velocity f , but this difference is more pronounced at large velocities. 

The relative location of the curves for given v and h and different a can be ex- 
plained by frame dragging. Namely, a string in the retrograde motion is effectively 
slowed down when it passes near the black hole and hence it spends more time near 
it. As the result its displacement is greater than the displacement for the prograde 
motion with the same impact parameter. Let us make order of magnitude estimation 
of this effect. For this purpose we use metric ( ^.1^ ). We focus our attention on the 
Lense-Thirring term and neglect for a moment the Newtonian part. Consider first a 
point particle moving in this metric with the velocity v and impact parameter h. In the 
absence of rotation it moves with constant velocity so that [X = vT, Y = b, Z = 0) 
and the proper time r is 

To = Vl-v^T . (5.1) 
In the presence of rotation one has 

dr^ = {l-v^)dT'-^dT\ (5.2) 
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Plot of K VS. impact parameter b for v/c=0.762 




Figure 9: Displacement parameter function of the impact parameter h for 

different velocities v. 
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where 

For large b one can write 
where 



R^ = b^ + v'^T^. (5.3) 
r = To + Ar , (5.4) 



_ 2Jbv r°° dT _ 4J 

(52 + ^2^2)3/2 - 5^/13^- ^ • ^ 

This quantity Ar characterize additional time delay for the motion in the Lense- 
Thirring field. The characteristic time of motion in the vicinity of a black hole is 
Tint ^ b/v. Thus we have 

. (5.6) 

One can expect that a similar delay takes place for a string motion. As a result of 
being longer close to the black hole the string has a larger displacement by the value 
Ak such that 

Ak Ar 4Ma . , ^ , ^, 

smh/3. (5.7) 



Tint b^ 

Quahtatively this relation explains dependence of k on a and /? presented in figure |. 

More generally, if we look at the a = a/M in the metric as a parameter, we can 
use perturbation theory to calculate the effect of rotation. The solution can be 
expanded in powers of a 

X'' ^aX^^o? , (5.8) 

i 

where X^{t, a) is the solution for the non-rotating case and X^{t, a) are the pertur- 
bation corrections. 

Numerical calculations show that the first two terms in the perturbation series 
give very good approximation to the solution as long as we are not close to critical 
scattering. As a consequence 

K = — Kia — , (5.9) 

where kq is the Y displacement for the non-rotating black hole, and /«i,/t2 are defined 
simply as 

fi:i,2 = lim A^'2(r,0) . (5.10) 

T— >00 

1,2 

In the above X'^ denotes the Y coordinate corrections. The equations of motion for 
1,2 1,2 

X^ are linear and X^ can be obtained in one run together with X^ ■ 

In order to determine the coefficients in equation ( |5.9|) we need to extrapolate 

the relevant data obtained during the simulation to infinity. One possibility is to 
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numerically advance the solution very far from the black hole and simply take the 
value at the end of the simulation as our estimate. The problem with this approach 
is that to advance very far by maintaining high accuracy of the solution is a very 
lengthy process, feasible only for larger impact parameters (since the grid can be 
sparser) . 

Our approach was to advance to a moderate distance (Xmax ~ 8000M) keeping 
high accuracy of the solution (judged by the constraint equations). Then we fit the 
calculated data by the function 



where k, k, Tgh, 7 are the parameters to be fitted. We take our first data point at 
X = lOOOM. To see how consistent this procedure is and to estimate the errors we 
perform three fits with different set of data points. The first set uses all the data 
points, the second one uses only the first half, and the third one uses only the last half 
of the data points. For this particular problem the disturbances from the boundary 
have adverse effect to the fit therefore we use the null boundary approach. 

Our experience shows that the fitting procedure for k works much better for larger 
velocities and larger impact parameters. For the velocities v/c = 0.762, 0.964, 0.995 
the estimated values of k from all three fits are consistent and differ by no more than 
3% 0. This fact gives us reasonable confidence that the plotted values of n indeed 
represent the true values. 

Unfortunately the situation is not so good for slower velocities. For example 
for the case v/c = 0.462 the slope of the fitted function k{0,t) is almost constant 
during the later phases of the simulation and fits using different sets of data points 
yield inconsistent (and improbable) values of k. We expect the behavior of k,{0, t) to 
change but we did not see it in our data even after driving the simulation four times 
further then for the other velocities (X^ax ~ 32, OOOM). Again, this behavior is more 
pronounced for smaller impact parameters h. 

To obtain the values for ni and K2 we use the same method as for k. We did not 
encounter any difficulties in estimating the values of /ti and /€2 for any of the four 
velocities. The trend here is the opposite of the one for k - the fits work best for 
smaller velocities and smaller impact parameters. The fits with different set of data 
points are typically less than 1% appart. 

Figure ^ shows the results in a log-log plot. Figure |TT] shows the plot of ki/k. 
The straight lines shown are obtained from least square fit. Note that the plots for 
the velocities v/c = 0.964 and v/c = 0.995 lie on top of each other. 

In the interval shown /ti and /€2 can be reasonably approximated by a functions 
of the form 



/t(0,r) 



7 



(5.11) 




(5.12) 



^For b > 12M the results from the three fits differ by less than 1%. 
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Plot of log,n(/f,/M) vs. log,n(b/M) 



Plot of log,(,(-ff./M) vs. log,„(b/M) 
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Figure 10: Plots of ki and k,2 for different velocities v and impact parameters b. The 
straight lines represent a linear fit. 

Plot of log.,g(/c.|//c) vs. log.|Q(b/M) 



1 ' r 

o v/c=0.995 

* v/c=0.964 

+ v/c=0.762 




1.2 1.4 
logio(b/M) 



1.6 



Figure 11: Plot of ki/k for different velocities v and impact parameters b. The 
straight lines represent a linear fit. 



Similarly for Ki/k, 



A3 



(5.13) 
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The values for 2,3 and Ai,2,3 obtained from a linear fit are shown in table |l[ By 
comparing ( |5.13| ) with (|5.7| ) one can conclude that the numerical value A3 = 2.23 is 
close to the value 2, which enters the relation ( ^.7|) . 

Note that the last data point on the plot for K2 is a bit out of a straight line. We 
are not sure if this is a genuine feature or simply inaccurate data points (e.g., due to 
round-off errors)^. We did not include these data points into calculation of the values 
m table I. 



v/c 


0.462 


0.762 


0.964 


0.995 


Ai 


22.3 


66.6 


231.0 


631.6 


Ai 


2.31 


2.36 


2.36 


2.35 


A2 


-83.9 


-381.4 


-1487.9 


-3717.0 


A2 


3.77 


3.89 


3.87 


3.83 


^3 




5.42 


6.21 


6.17 


A3 




2.23 


2.23 


2.23 



Table 1: Values of fitted parameters from equation ( ^.12| ) and (|5.13|) 



5.2 Form of the kinks 



Profile of the kink for b=80M,v/c=0.762,a/M=0 



Profiles of the kinks for b=6M,v/c=0.762,a/M = - 1 ,0, 1 
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Figure 12: Profiles of the kinks for different impact parameters 



Figure 12 shows details of the transition from the "old phase" to the "new phase" 
for two different impact parameters. We see that the width of the kinks differs 



^The discrepancies are larger than those which could result from the extrapolation procedure. 
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significantly. In general the width also depends on the velocity v/c The dependence 
of the kink width on velocity and the impact parameter can be estimated from the 
weak field approximation as 

he 

w ~ — . (5.14) 

V 

We operationally define the width of the kink to be the width of the peak of dY/ da 



at l/20th of its hight. Figure [T^ shows comparison of the kink widths obtained from 
simulation with a function 

he 

Wanalyt = ''•l^— ■ (5.15) 

The numerical constant was obtained by a fit. The match with the data is very good. 

We also noticed that the dependency of the width of the kink on the rotation 
parameter a is rather small. In general 



a/M=l < Wa/M=0 < Wa/M=-1 ■ {5.16) 



For the parameters shown on figure |T3| the rotation made a difference from 3M to 
7M for the extremal black holes. 



6 Discussion 

In this paper we studied the scattering of a long test cosmic string by a rotating 
black hole. We demonstrated that qualitatively many general features of the weak 
field scattering are present in scattering of the cosmic string in the strong field regime. 
Displacement of the string in the y— direction always has the form of a transition 
of the string from the initial phase (initial plane) to the final phase (final plane 
which is parallel to the initial one and displaced in the direction to the black hole by 
distance k). The boundary between these phases is a kink moving with the velocity 
of light away from the center. An important difference between weak- and strong- 
field scattering is in the dependence of k on the impact parameter h and velocity v. 
In general k for strong field scattering is much greater that its value calculated by 
weak-field approximation. It is also always greater for retrograde scattering than for 
a prograde scattering. The explanation of this is quite simple: The retrograde string 
spends more time in the vicinity of the black hole than a prograde one. This effect is 
a result of the dragging of the string into rotation by the black hole. 

The width of the kink for the strong field scattering is smaller than for the scat- 
tering in the weak field regime. 

As was demonstrated in previous studies for string scattering by non-rotating 
black holes, for any given velocity there exists a critical value hcrit of the initial impact 
parameter so that for h < hcrit the string is trapped by the black hole, while for h > hcrit 
it is scattered by the black hole (and all points on the string reach infinity). In this 
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Figure 13: Comparison of kink widths for v/c = 0.462,0.762,0.964,0.995 and im- 
pact parameters h/M = 9, 12, 15, 18, 25, 40 obtained from simulation (circles) and the 
analytic approximation ( |5.15|) . The data shown are taken for a/M = 0. 



paper we focus our attention on scattering in the regime where b is significantly greater 
than bcrit- The study of string capture and near-critical string scattering requires 
modification of the computational program and more time consuming calculations. 
We are going to present these results in a subsequent publication. 
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A Numerical scheme details 



In section |3.2| we briefly described the numerical scheme we adopted. In this appendix 
we present some more details. Figure ^ shows the overall structure of the time domain. 

The numerical grid in the sigma direction is non-uniform. This is very important 
since there are regions where the string is relatively straight and regions with high 
curvature. By far the numerically most "sensitive" region is the kink. Therefore we 
make the grid denser in the vicinity of the kink. The length of the dense zone is 
chosen to be approximately twice the kink widthQ. The steeper the kink the denser 
the grid must be. Typically, the dense zone is 3-30 times denser then the rest of the 
grid. 

The time step At must not exceed the smallest distance between any two grid 
points Act. A convenient choice turns out to be 



At = Aay/l-vyd^ . (A.l) 

This means that dense grid also implies small time step which in turn makes the 
simulation time longer. The formula (|5.14| ) explains why it is more difficult to deal 
with small impact parameters and relativistic velocities. 

Since the kink is moving away from the central point a = the dense zone must 
move with it. Therefore we regularly check the position of the kink and adjust the grid 
when the kink reaches the boundary of the dense zone. This procedure is illustrated 



in figure Note that we introduced similar (but much shorter) dense zone at the 
edge of the string. This is to prevent from "cutting" the string too quickly (every 
time step we lose one grid point). When we change the the grid structure we must 
use interpolation to obtain the variable values at the new grid points. In our case we 
used cubic spline which proved to work very well. 

At certain stage the two dense zones merge and eventually they disappear alto- 
gether (which means that the kink passed the edge of our numerical domain) . At this 
point we can increase At to refiect the sparser grid. 

Since our problem is naturally formulated on rectangular grid the finite difference 
method is used to discretize the equations of motion. To approximate the first and 
second derivatives we use standard second order centered formulas, e.g., 

{dlxn., = ^"^'^ , (A.3) 

where index i marks the z-th time step and j marks the j-th grid point in the sigma 
direction. Similar expressions are used for the spatial derivatives except for the points 



*The grid within each zone is uniform 
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END OF SIMIILATION 




^ » DENSE ZONE f 

G 0=0 

Figure 14: Structure of the numerical grid 

lying at the boundary between dense and "normal" zones. In that case we must use 
modified formulas 

(A.5) 

where A(Tj_i = aj — aj-i and Actj+i = Uj+i — Uj. 

After the discretization process we obtain a set of four coupled non-linear equa- 
tions for the four unknown variables at each grid point aj. Since the equations 
for different j's are independent of each other we can solve them in parallel. 
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